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ABSTRACT 

We report measurements of the carbon monoxide ground state rotational transition ( 12 C 16 J = 1—0) 
with the Zpectrometer ultra- wideband spectrometer on the 100 m diameter Green Bank Telescope. 
The sample comprises 11 galaxies with redshifts between z = 2.1 and 3.5 from a total sample of 
24 targets identified by Hersche '-ATLAS photometric colors from the SPIRE instrument. Nine of 
the CO measurements are new redshift determinations, substantially adding to the number of de- 
tections of galaxies with rest-frame peak submillimeter emission near 100 fim. The CO detections 
confirm the existence of massive gas reservoirs within these luminous dusty star-forming galaxies 
(DSFGs). The CO redshift distribution of the 350 /mi-selected galaxies is strikingly similar to the 
optical redshifts of 850 /zm-selected submillimeter galaxies (SMGs) in 2.1 < z < 3.5. Spectroscopic 
redshifts break a temperature-redshift degeneracy; optically thin dust models fit to the far-infrared 
photometry indicate characteristic dust temperatures near 34 K for most of the galaxies we detect in 
CO. Detections of two warmer galaxies and statistically significant nondetections hint at warmer or 
molecule-poor DSFGs with redshifts difficult determine from from Herschel-SPIRE photometric colors 
alone. Many of the galaxies identified by H-ATLAS photometry are expected to be amplified by fore- 
ground gravitational lenses. Analysis of CO linewidths and luminosities provides a method for finding 
approximate gravitational lens magnifications /i from spectroscopic data alone, yielding /i ~ 3-20. 
Corrected for magnification, most galaxy luminosities are consistent with an ultra-luminous infrared 
galaxy (ULIRG) classification, but three are candidate hyper-LIRGs with luminosities greater than 
10 13 L Q . 

Subject headings: galaxies: high redshift — galaxies: ISM — galaxies: evolution — submillimeter: 
galaxies 
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1. INTRODUCTION 

Observations of the far-IR/submillimeter background 
with CO BE demonstrated that a substantial fraction of 
the uni verse's star formation took place beh ind a veil 
of dust (jPueet et al.lll996t iFixsen et al.lll998[ ). Because 
these were integrated measurements, however, they could 
not identify which populations of dust-obscured galaxies 
contained this vigorous star formation. A breakthrough 
in resolving the background came in the late 1990s, with 
imaging by the Submillimctcr Common-User Bolometer 
Array (SCUBA) on the James Clerk Maxwell Telescope 
(JCMT). Its initial surveys at 850 urn (ISmail et al.l ll99ft 
iBarger et all 119981 : lHugh.es et all ll99S|) detected a new 
population of bright (> 5mJy) galaxies. Named af- 
ter the wavelengths where they are most visible, these 
submillimeter galaxies (SMGs) are systems with appar- 
ently vast (> 10 13 L©) bolometric luminosities but with 
such high obscurations that their optic al counterparts 
are faint or absent (see lBlain et al.ir2002L and references 
therein). SMGs brighter than SCUBA's confusion limits 
could not account for all of the 850 /im background, but 
clearly made a substantial contribution to it. 

Over the last fifteen years, much of the effort to un- 
derstand the origins of the far-IR/submillimeter back- 
ground has focused on bright SMGs selected from 850- 
1200 /im surveys. SMGs are sometimes treated as rep- 
resentatives of a more general category of high-redshift 
galaxies, dusty star-forming galaxies (DSFGs) whose lu- 
minosities are dominated by obscured star formation. A 
major initial hurdle was verifying that bright SMGs actu- 
ally do lie at high redshifts. W hile two early SCUBA de - 
tections had optical redshifts ([Ivison et al.l[i"9 98. 2000), 



quick l y con firmed with CO spectroscopy (jFraver et al.l 
[19981 Tl999h . progress in measuring redshifts of other 
SMGs fou ndered due to their v ery high obscurations. 
Only after iChapman et al.l (|2003l ) took advantage of ra- 
dio continuum mapping to determine precise positions 
for blind optical spectroscopy did it be come possible to 
obtain CO detect i ons in large number s ([Neri et alj|2003t 
iGreve et al.l 120051 : iTacconi et al.ll2006ft . and to derive an 
SMG redshift distribu tion peaking around z ~ 2-2.5 
(jChapman et al.l [2005) . Heroic efforts to explore the 
high-z tail that radio pre-selection misses when radio 
counterparts fall below typical survey flu x limits have 
ident i fied a handful of sour c es at z > 4 (e.g.jCapak et all 
2001: ICoppin et all l200l I2010bt IRiechers et al.l 120101: 



Capak et al.l 120 111 ), but detailed analysis limits the pos- 
sible significance of this tail in the ~ 1 mm population 
([Ivison et alJl2005t IWardlow et al.H201lh. We now know 



that bright SMGs h ave large stellar (ISmail et al 



Hainline et all 12011ft . molecular gas (IGreve et al 
Tacconi et al.ll2008ft . and dynamical (iGenzel et al 
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mass es, that many of them are mergers ([Conselice et all 
2003), and that their large l uminosities are powered prin- 
cipally by star fo rmation ([Alexander et al.l 120031 120051 : 
iPope etall 12008ft . Explaining the observed properties 
of bright SMGs, their evolutionary states, and their re- 
lationships to populations of galaxies selected at other 
wavelengths is a cur rent major challenge for galaxy evo- 
lution models (e . g.. iBaueh et al.l 120051: i Swinbank et all 
l200l lDavtetafll201tt: iSomerville etallboill ). 

In parallel with the growth in our understanding of 
bright SMGs, it is also becoming clear that current sam- 



ples give an incomplete picture of the full variety of 
DSFGs. First, 850 /im sources fainter than SCUBA's 
nominal confusion limit of about 2 mJy, although ac- 
cessible via gravitational lensing (e.g.,|Smail et alJl2002t 
iKneib et al.lT2005t iKnudsen et all 120081 must at some 
point start to resemble optically sel ected galaxies more 
than heavily obscured SMGs (e.g., iGenzel et all 120101 : 
IDaddi et all l20loh . Second, even among bright DS- 
FGs, 850 /im-bright SMGs have distinct selection biases. 
DSFG samples selected at longer wavelengths appear 
to have cooler dust temperatures and higher median 
redshifts (iDannerbauer et all 120041: j Valiantc et all l~ 



IWardlow et al.l 1201 H ILindner et al.ll20Ti!T DSFG sam- 
ples selected at shorter wavelengths, conversely, include 
populations wit h warmer dust that 850 /im-selected sam - 
ples can miss, (Blain et al. 2004; IChapman et al.l[2004ft . 
and that tend to have both lower redshifts and more 
bolo metrically significant active galacti c nucl ei (AGNs; 
e.g-JHouck et all 120051: iWeedman et al.ll2006UYan et all 
12007ft . Third, existing samples of SMGs with spectro- 
scopic redshifts often suffer from biases associated with 
the steps used to determine those redshifts. For exam- 
ple, precise localization of a submillimeter source usu- 
ally relies on radio continuum mapping, and while the 
deepest VIA imagi ng yields counterpa rts for a high 
fraction of DSFGs (jLindner et alll201lft . more typical 
VLA maps tend to deliver counterparts for only 60- 
70% of SMGs. Subsequent optical spectroscopy based 
on these positions fail s to yield redshifts fo r a modest 
fraction of candidates ([Chapman et al.l 1200 5). and even 
when apparently successful, attempts to obtain CO de- 
tections of the gas reservoirs associated with these mas- 
sive starbursts can sometimes fail to yield conclusive 
confirmation, r aising questions oy er either the identifica- 
tion or redshift ([Greve et al.ll2005ft . Mid-I R spectroscopy 
can avoid some of these difficult ies (e . g.. |Valiante et all 
200ft iMenendez-Delmestre et all[2007l: IPope et alll200S 



Menendez-Delmestre et all l200a ICoppin et al.l 12010a 
but suffers from its own problems for sources with power- 
law spectra or confusion from multiple sources within 
large beams or slits. Finally, many of the seminal stud- 
ies of SMGs mapped relatively small areas on the sky. 
Notwithstanding the large line-of-sight interval probed 
by 850 /im selection, which can partly compensate for a 
small area, SMGs are such rare sources (mergers caught 
at special moments, with the most luminous caught at 
the most special of moments) that cosmic variance re- 
mains a concern for the derived redshift distributions. 

New instruments capable of producing deep images of 
large regions of the sky, and of conducting efficient spec- 
tral surveys over wide band widths, have accelerated the 
discovery of high-redshift DSFGs with a broader range 
of physical properties than could be probed by previ- 
ous efforts. Surve y areas at A ~ 1 mm have increased 
substantially fe.g.lMarsden et al|[2009HWeifi et al.l l2009t 
lAustermann et aill2010l ; lHatsukade et al.ll2011ft . with ex- 
tremely wide-area surveys possible both from the ground 
([Vieira et alll2010l : lMarriage et al.ll2011l ). and from space 
with the Herschel Space Observator^l (Pilbratt et al. 
2010). Large-area surveys are identifying many very 

24 Herschel is an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator consortia 
and with important participation from NASA. 
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bright DSFGs whose fluxes are rivaled by those of only 
a few extreme, serendipitously discovered objects that 
have been confirmed to be gravitationally lensed (e.g., 
ISwinbank et al.l l2010T l . Regardless of the balance be- 
tween intrinsically hyperluminous systems and less ex- 
treme but gravitationally lensed galaxies within these 
samples, we are no longer missing the rarest DSFGs be- 
cause of limited sky coverage. Herschel is playing a par- 
ticularly important role because the wavelen gth coverage 
of its SPIRE instrument ([Griffin et al.ll2010D allows selec- 
tion of DSFG sa mples that are relativ ely free of dust tem- 
perature biases ([Magdis et al.ll20Tol ). although they are 
limited by confusio n to only the most extre me luminosity 
systems at z > 1 (|Svmeonidis et al.ll2011l) . Specialized 
instruments designed for wide-band spectral line surveys 
now enable the determ i nation of blind C O redshifts (see, 
e.g.. IBaker etHI [20071: IWeifi et al.ll2009h for bright DS- 
FGs, without intermediate radio continuum mapping or 
optical spectroscopy. An example of the combination 
of these new developments is the recent use of two wide- 
bandwidth instruments to obtain CO redshift s for a com- 
plete sample of five brigh t Herschel sources ([Lupu et al.l 
120101: IFraver et al.l 120111 ) , an essential st ep in confirm- 
ing t heir status as galaxy-galaxy lenses ([Negrello et al.l 
120101) . 

Expanding on this initial work, we here report on A ~ 
1 cm spectroscopy of the 12 C 16 J = 1 — ground-state 
rotational transition toward two dozen of the bright- 
est DSFGs in catalogs from the Herschel A strophysical 
Terahertz Large Area Survey (H- ATLAS; lEales et al.| 
201( 1 llbar et al.1 l20Tot iPascale et al.l I20111 iRigbv et al.l 
201ltlSmith et al.ll20lHT program, using the Zpectrome- 



ter ultra-wideband spectrometer on the National Radio 
Astronomy Observatory's 100-meter diameter Robert C. 
Byrd Green Bank Telescope (GBT). Submillimeter con- 
tinuum flux ratios provide some coarse redshift informa- 
tion for many sources, but the precise redshifts needed to 
enable additional observations with narrow-band instru- 
ments require spectroscopy of atomic or molecular lines. 
The Zpectrometer is one of several ultra-wideband spec- 
trometers built for this purpose, and is the first instru- 
ment to make routine measurements of the CO J = 1 — 
rotational transition from high-redshift galaxies (e.g. , 
Swinbank et~all 120101; lHarris et al.l [20101; [Negrello et al 



2010t[Fraver et al.ll201 ll: iScott et alil201lURiechers et al 
2011cl lbD With the CO molecule's J = 1 level only 5.4 K 



above the ground state, its low but nonzero permanent 
dipole moment, and its strong C-0 bond, the ,7=1 — 
transition is the best tracer of molecular gas over a wide 
range of conditions in molecular clouds. In addition to 
giving a spectroscopic marker for redshift measurements, 
velocity-resolved spectroscopy yields dynamical informa- 
tion, which together with gas masses derived from CO 
intensities provides key inputs to understanding DSFGs' 
star formation efficiencies and overall evolutionary states. 

Subsequent sections of this paper describe the obser- 
vations and the initial results from this sample. Sec- 
tion [2] describes target selection and observations, Sec- 
tion [3] contains observational results, and Section [4] pro- 
vides further analysis and discussion. Calculations use 
a ACDM cosmology with = 0.27, Q\ = 0.73, and 
h = 0.71 (ISpergel et al.ll2007l) . 



2. OBSERVATIONS 

We drew our targets from early Herschel- ATLAS cata- 
logs (H-ATLAS collaboration, priv. comm. 2010) of con- 
tinuum detections in the Herschel SPIRE instrument's 
250, 350, and 500 /jm wavelength photometric channels. 
For this initial study, we selected bright "350 /im peaker" 
galaxies with flux densities 5^(350 /zm) > 115 mJy and 
observed spectral energy distributions (SEDs) peaking 
in the SPIRE 350 /im band, within errors. The cata- 
logs screened out local spiral g alaxies and high-z bl azars 
with the methods de scribed in lNegrello et al l (|2007f ) and 
iNegrello et al.l ([20101 ) . This provided a set of bright tar- 
gets with peak far-infrared emission near 100 /jm wave- 
length in the rest frame for galaxies with z ~ 2 to 
3. Most of the targets were in the equatorial multi- 
wavelength Galaxy And Mass Assembly (GAM A) survey 
(|Driver et al. 2009) fields near Right Ascension 9, 12, and 
15 hours, with a few targets from the North Galactic Pole 
(NGP) field near a = 13 hours and 5 = 27°. Table [I] is a 
list of target positions and integration time information. 

We observed the targets with the Zpectrometer ana- 
log lag cross-correlator spectrometer ([Harris et al.l 120071 : 
lHarrisI 120051) connected to the GBT's facility Ka-band 
receiver, which was configured as a correlation re- 
ceiver. Receiver improvements in Fall 2010 extended its 
high frequency performance, allowing spectroscopy from 
25.6 to 37.7 GHz, corresponding to the CO molecule's 
115.27 GHz J = 1 - transition at redshifts of 2.1 < 
z < 3.5. The Zpectrometer's spectral resolution is a sine 
function with full-width at half maximum (FWHM) of 
20 MHz, sufficient to provide a few resolution elements 
across typical galaxy lineshapes. Given the instrument's 
wide overall bandwidth, the spectral resolution varied 
with frequency across the spectra from 234 to 157 km s _1 , 
and the FWHM beamsize from 27 to 16 arcsec. 

A combination of the correlation receiver's ability to 
difference power between two beams on the sky and the 
Zpectrometer's large bandwidth allows different observ- 
ing techniques from those common in narrowband total- 
power radio astronomy. Our instrument, observing tech- 
nique, data reduction, and calib ration methods are fully 
described in lHarris et al.l (|2010[ ). Briefly, the correlation 
receiver implementation electrically differences power be- 
tween the receiver's two input horns, which are sepa- 
rated by 78 arcsec on the sky. We switched the source 
between the two horns by moving the GBT's secondary 
mirror 78 arcsec in a 10 s cycle, then differencing spectra 
from the two positions to obtain a source spectrum. To 
eliminate the residual tens of mJy of spectral baseline 
structure from optical offsets, we observed pairs of tar- 
gets close in position on the sky, cycling between sources 
every 4min, again differencing this pair. In this differ- 
ence spectrum of the two positions an emission line from 
the first source would appear in the positive sense, while 
one from the second source would appear in the negative 
sense. Differencing left little optical offset and baseline 
structure, at the cost of eliminating information about 
individual source continua. 

Observations were conducted in sessions of 3 to 8 hours 
duration (see Table [T]) on dates from 2010 November 
through 2011 April for a total of 64.5 hr of on-sky ob- 
serving time. Counting observing overheads, the actual 
elapsed observing time was about 130 hr. We observed 



4 



Table 1 

Target list grouped by observed pairs (see Sec.[2]l. 
Positions are encoded in the source names following the 
IAU convention of providing Right Ascension and 
Declination, here in J2000.0 coordinates. Other 
columns give the GAMA or North Galactic Pole (NGP) 
field, total integration time, number of individual 
observing sessions, and a target reference number for 
figures and other tables; letters denote targets with 
redshifts f rom the H- ATLAS Science Demonstration 

Phase UNegrello et aTll201CI ; a is ID. 17b, see also 
ILupu et al,ll201Gt b is ID. 130. see also lFraver et ahl 
120111) . Emission from the first source in each pair would 
appear in the positive sense in the spectra in Figure [T] 
and emission from the second source in the negative 
sense. Positions and therefore names of sources in the 
NGP field are preliminary and may change slightly. 



H-ATLAS 


Field 


tint 


No. Target 






[hr] 


sess. No. 



1 

2 



3 
4 

a 
b 

5 
6 

7 
8 

9 
10 

11 
12 

13 
14 

15 
16 

17 
18 

19 
20 

21 
22 



in a variety of weather conditions, mostly with reason- 
able to good Ka-band atmospheric transmission and low 
wind. 

We established absolute flux scales across the spectra 
by dividing the astronomical source difference spectra 
by the spectrum of a bright (few Jy) quasar suitable 
as a pointing reference near each target pair. In some 
cases we could use one of the cm-wave flux standards 
3C48, 3C286, or 3C 147 directly (0.80, 1.83, 1.47 Jy at 
32 GHz, respectively; iThe Astronomical Almanac] 2011), 
but for most of our sources we cross-calibrated spec- 
tra of the nearby pointing source with one or more 
of the s tandards. We cross-checked 3C48's flux den- 
sity from lThe Astronomica l Almanad against a Mars flux 
density model bv B. ButleJ^I and found agreement within 
2%. Comparison with the recommended 2012 January 
Jansky-VLA flux density for 3C286 is 1.96 Jy (7% higher) 

25 http:/ /www. aoc.nrao.edu/~bbutler/work/mars/model/ 



than the value we use, and we find a Ka-band spectral in- 
dex of a = —0.8 instead of the Jansky-VLA's a = —0.4. 
Cross-calibrations of the flux density for the quasar we 
used for the GAMA09 sources against all three of the 
standards agreed within 10%. We determine calibrator 
spectral indices across the Ka-band, which can be impor- 
tant in transferring the 32 GHz standard fluxes to other 
frequencies in the band, by comparison with Mars' black- 
body a = +2. 

We pointed and focused every hour on the bright 
quasars we had selected as secondary flux calibrators 
near our sources. Pointing offsets were always within a 
third of a beam near band center. We took spectra of the 
pointing source at the beginning and end of each pointing 
cycle to measure changes in optical gain and atmospheric 
transmission. Dividing the difference spectrum of an as- 
tronomical source pair by the average spectrum of the 
pointing source not only corrected for bandpass gain but 
also for the effects of pointing errors and changing atmo- 
spheric transmission with timescales comparable to an 
hour, or longer. Session to session repeatability for the 
quasar flux densities was generally within 10%, and we 
take 20% as representative of the overall amplitude cali- 
bration uncertainty to include gain effects from pointing 
drift. 

Version D of the Zpectrometer's standard GBTIDL 
scripts provided quick-look assessment during observa- 
tions and produced files for combining data from different 
sessions. We used version 5. 4 of our Zred package, which 
is wr itten in the R language (|R Development Core Teaml 
2006), for further analysis. The receiver produced sys- 
tematic spectral structure with fluctuating amplitude 
across many of the spectra. Spectral lines from galax- 
ies were much narrower than features in the systematic 
structures, which enabled us t o remove the instrumen- 
tal artifacts ([Harris et al.ll2010l contains further informa- 
tion). We removed a systematic ripple with Fourier fil- 
tering and a complex but fixed-pattern structure by fit- 
ting and subtracting median-filtered session-average tem- 
plates of the structure. We removed no further base- 
line structure, and with these corrections we could keep 
nearly all of the data. 

2.1. Line detection algorithm 

In this exploratory phase of a full survey, we set inte- 
gration times to identify bright CO sources rather than 
to get high signal-to-noise ratios on individual sources. 
Since system noise changes across the Zpectrometer's 
12.1 GHz band and nonideal noise is clearly present, at- 
tempting to estimate noise for signal-to-noise ratios by 
computing the channel fluctuations across the spectrum's 
entire frequency range would overestimate noise in some 
parts of the spectrum and underestimate it in others. 
Instead, we harness time-series information from sub- 
integrations to make individual channel no ise estimates 
as par t of a detection confidence algorithm. lHarris et al.l 
(2010) has a complete discussion of this algorithm, which 
analyzes whether the amplitude in each frequency chan- 
nel is statistically higher than the average of its neigh- 
bors, based on a generalized Student-t test of many boot- 
strapped realizations of each spectrum. 

The algorithm's main advantage over traditional meth- 
ods is that it treats noise in individual frequency chan- 
nels. Its main weakness stems from its assumption that 



J083051. 
J084933. 

J083929. 
J084259. 

J090302. 
J091305. 

J091840. 
J085111. 

J091948. 
J092135. 

J113526. 
J113243. 

J113833. 
J113803. 

J114637. 
J115112. 

J115820. 
J114752. 

J132426. 
J133008. 

J134429. 
J133649. 

J141351. 
J142751. 



0+013224 
4+021443 

5+023536 
9+024958 

9-014127 

0- 005343 

8+023047 
7+004933 

8- 005036 
6+000131 

3-014605 

1- 005108 

3+004909 
5-011735 

9- 001132 
3-012638 

2- 013753 
7-005832 

9+284452 
3+245860 

4+303036 
9+291801 

9-000026 
0+004233 



GAMA09 
GAMA09 
GAMA09 
GAMA09 
GAMA09 
GAMA12 
GAMA12 
GAMA12 
GAMA12 
NGP 
NGP 
GAMA15 



7.87 
7.08 
5.46 
4.59 
3.67 
5.64 
3.67 
7.35 
5.25 
2.89 
3.94 
10.23 



5 



the underlying spectral baselines are relatively smooth 
on small scales so that local averages are representative. 
Failure of this assumption keeps the computed confidence 
level from being absolute, as discontinuities in spectral 
baselines and small fluctuations near bright features can 
all register as potential detections. 

Our experience with the algorithm on Zpectrometer 
data is very good, based on comparisons of tentative de- 
tections within subsets of data and from detections of 
different lines from the same galaxies with other instru- 
ments. Given deviations from the assumptions caused 
by nonideal noise of a few hundred /iJy, the algorithm 
remains a powerful guide rather than an exact indica- 
tor for identifying weak lines. Experience gained from 
viewing many spectra toward different sources under dif- 
ferent weather conditions, but all with the same fre- 
quency coverage, provides valuable impressions of typical 
baseline structures and noise behavior, further aiding in 
screening against spurious detections. All detections here 
come from at least two independent observing sessions, 
and most have at least tentative detections in individ- 
ual sessions. By generating 1000 bootstrapped spectra 
as part of the detection algorithm, the data are thor- 
oughly mixed by time and session. Sweeping through a 
wide range of binnings ensures that a detection is not 
based on a single favorable choice. We have been con- 
servative in our decisions, while recognizing that a small 
number of false positives are more beneficial than false 
negatives in observations designed to spur followup work. 
In practice, independent observations of other lines from 
the same targets have justified our detection selections. 

3. RESULTS 

Panels in Figure [T] show the spectra and confidence 
plots for all of the source pairs we observed. The up- 
per panel for each source pair is the spectrum across the 
Zpectrometer band. Vertical dashed lines mark line fre- 
quencies of detected galaxies. The lower panel shows the 
detection probabilities p from our detection algorithm, 
versus frequency, given as Confidence — — log 10 (l — p) 
(numerically, the scale is equivalent to the number of 
nines in confidence: 0.9, 0.99, 0.999, etc. for p > 0.9). 
Each dot in the plot is an individual channel's confi- 
dence measure (within the algorithm's assumptions) for 
a given combination of binning width and starting point. 
Columns of dots at specific frequencies show where a po- 
tential detection is relatively immune to exact binning 
parameters, indicating that a real line is present rather 
than a favorable binning for a chance fluctuation. 

Our brief comments on individual spectra are: 
J083051. 0+013224 and 30849334+021443 (Fig.[TJ 1-2): 
Only the second target in the pair (emission appears 
in the negative sense) is detected. D.A. Riechers (priv. 
comm. 2011) has detected a single strong line with the 
CARMA observatory toward the first in the pair. If the 
line were CO J = 3 - 2, J083051.0+013224 would lie 
at a redshift with the 1 — line in a low-noise region of 
the Zpectrometer band. This is a statistically significant 
nondetection that we discuss in section l3~Tl 
J083929. 5+023536 and J 084259. 9+ 024958 (Fig. DO 3-4): 
Neither target is clearly detected in this spectrum. 
3090302.9-014121 and 3091305.0-005343 (Fig.EJ a-c): 
Clear detections of both targets. The high confidence 
at slightly lower frequency than the strong positive line 



is likely an artifact, as even a modest dip can be far 
from the local amplitude mean, a signature the detection 
algorithm interprets as a line. 

3091840.8+023047 and 3085111.1+004933 (Fig. ED 5-6): 
Clear detection of the first target in the pair. The con- 
tinuum offset from zero shows that the first target has a 
higher continuum flux than the second target. 

3091948.8- 005036 and 3092135.6+000131 (Fig. [Q 7- 
8): Neither target is clearly detected in this spectrum. 
Residual large-scale structure may obscure what could 
be tentative detections. Even small noise fluctuations at 
the tops of large-scale positive and negative structures 
are far from local means, so they register strongly in the 
confidence plot. 

3113526.3-014605 and 3113243.1-005108 (Fig. [I] 9- 
10): Clear detections of both targets. Some spurious 
high confidence peaks are associated with each of the 
bright lines. 

3113833.3+004909 and 3113803.5- 011135 (Fig. [TJ 11- 
12): Neither target is clearly detected in this spectrum. 

3114631.9- 001132 and 3115112.3-012638 (Fig. [U 13- 
14): Clear detection of the first target. 

3115820.2- 013153 and 3114152.1-005832 (Fig. [U 15- 
16): Clear detection of the first target. 
3132426.9+284452 and 3133008.3+245860 (Fig. [Q 17- 
18): Tentative detection of the first target. Baseline 
structure to slightly higher frequencies makes it difficult 
to find a local baseline, so the line parameters are uncer- 
tain. 

3134429.4+303036 and 3133649.9+291801 (Fig. Q] 19- 
20): Clear detections of both targets, with one line in 
the positive sense and two in the negative sense. We at- 
tribute the stronger negative line with the SPIRE source. 
For completeness, we assign the other negative line a ten- 
tative formal detection because it is uncharacteristically 
broad for a spurious signal, but it is close to the positive- 
sense line where the local mean changes rapidly and may 
throw off the detection algorithm. 

3141351.9-000026 and 3142151.0+004233 (Fig. [U 21- 
22): A strong detection of the first target. The high 
confidence measures for dips to either side are most likely 
due to a high local mean in the region. The dips are 
rather wide to be due to a galaxy, so it is unlikely that 
they represent detections. The continuum slope indicates 
a difference in spectral index between the two targets. 

Overall, we detected 11 of the 24 targets in our 
sample. Two of the detections were blind indepen- 
dent confir mations of source s with establis hed redshifts 
(ID. 1 7b in lLupuet aT1[20Tol and ID. 130 in iFraver et all 
1201 If ). This success rate is similar to that of detec- 
tions from the Plateau de Bure millimeter-wave inter- 
ferometer starting from optical redshift catalogs (e.g. , 
Neri et all [200l iGreve et~aTI [20051 iTacconi et atl [20061 : 
Bothwell et al.ll2012h . but without the complications as- 
sociated with fin ding optical redshifts for submillimeter 
sources (see, e.g.. lChapman et al"1l2005l ). Section con- 
tains a more extensive discussion of detection complete- 
ness. In addition to the 11 detections, we also list two 
tentative detections, denoted by italics in Table [3] and 
(for J132426.9+284452) with open circles in the figures. 
Line parameters for tentative detections were too uncer- 
tain for robust error estimates. 

Based on the submillimeter photometric selection and 
line strength, we initially derived redshifts assuming that 
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Figure 1. The first twelve sets of spectra and detection confidence plots for targets labeled by target number in Table ITl and described in 
Section [3] The upper panel in each set is the difference spectrum for each pair (panel letter in the upper right hand corner). The lower 
panel is a detection confidence plot as described in the text. Coherent columns of dots indicate detections that are insensitive to binning 
parameters. Vertical dashed lines mark detections or tentative detections. 



the lines were the redshifted CO J = 1 — transition 
rather than CO J = 2 — 1 from a 5.1 < z < 8.0 galaxy 
or a line from a species other than CO. This assumption 
has proved correct for all 11 sources, which have other ob- 
served lines, most starting with redshifts from the Zpec- 
trometer observations (D.A. Riechers, priv. comm. 201 1; 
P.P. van der Werf, priv. comm. 2011: ILudu eT^MHl . 

Table [3] summarizes observed and derived source pa- 
rameters. For detected lines, the parameters are from 
single-component fits to Gaussian lineshapes, with er- 
rors given by the statistical uncertainties in the fit at 
the 68% ("lrj") confidence level. Gaussian fits to the 
convolutions of Gaussian line shapes and the correlator's 
sine instrumental profile shows that linewidth corrections 
are unimportant for linewidths above 200kms _1 FWHM 
(the correction is 12% when the Gaussian line and sine 
FWHMs are equal, falling below 1% when the Gaussian's 
FWHM is 1.5 times the sine's FWHM or wider). All 
lines are broader than this, so we make no corrections. 
The table also contains estimates of the total infrared 
(8-1000 /im) fluxes and dust temperatures obtained from 
fits to the SPIRE photometry (H- ATLAS collaboration, 
priv. comm. 2010), as discussed later. In addition to the 
detections from this program, Table [3] also includes CO 
1 — data from a previous Zpectrometer detec tion of H- 
ATLAS J090311. 6+003906 dFraver et alll27nil ID.81). 



3.1. Flux relationships and detection completeness 

We have firm detections of 11 of the 24 targets in our 
program, a number large enough to draw some sample 
conclusions. Here we use the molecular and continuum 
flux information to explore the detection completeness 
and to evaluate possible reasons for CO nondetections. 
Figure explores the relationship between 350 /im flux 
density, 5^(350), and the CO J = 1 — integrated flux, 
5 !y (CO)Aw. Sources without CO detections have zero 
amplitude in this plot, and the open circle denotes a 
tentative CO detection. Line nondetections generally 
fall to lower 350 /im flux densities, but several galaxies 
are bright in continuum but not detected in CO. We 
can take the third-brightest 5^(350) source in our sam- 
ple (J083051.0+013224, target number 1 in figures and 
tables) as an example of a nondctcction that implies a 
galaxy with a redshift outside the Zpectrometer band 
or an abnormally low CO to continuum flux ratio. The 
former seems more likely in this case: as noted above, 
CARMA has a clear detection of a mid- J CO line from 
this source, but deep integrations by the Zpectrometer 
and other instruments have not found other lines corre- 
sponding to a redshift within the Zpectrometer's z = 2- 
3.5 band. 

Most likely, some fraction of the nondetected galax- 
ies have redshifts that are not well predicted by con- 
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Figure 1 (cont.). The second twelve sets of spectra and detection confidence plots for targets labeled by target number in Table[T]and 
described in Section [3] 



tiiTuum properties, while others have CO 1 — lines 
that are fainter than our detection threshold. Before 
we simply ascribed CO nondetections to galaxies with 
observed luminosities below some threshold, we first ex- 
amined the line detection completeness and the rela- 
tionship between line and continuum fluxes. Estimat- 
ing line detection completeness in broadband spectra is 
complicated by system temperature and nonideal noise 
that vary with frequency; this frustrates any attempt 
to define a quantity such as the baseline rms in nar- 
rowband spectra that could specify a simple detection 
limit. We estimated completeness levels for nondetec- 
tions by simulation, adding sets of synthetic lines as 
frequency combs across subscans for all sources, run- 
ning the modified spectra through the data reduction 
pipeline, and inspecting the spectra to see which syn- 
thetic lines we could clearly identify. With a comb of 
seven 400 km s -1 wide lines (a width close to the me- 
dian of the astronomical source lincwidths) across each 
spectrum, we recovered 90% of the synthetic lines with 
CO integrated intensity S„(CO)Av = 800 mJykms -1 , 
70% of the lines with S„(CO)Av = 600 mJykms" 1 , and 
35% of those with S v (CO)Av = 400 mJykms- 1 . Taking 
S„(CO) Av — 600 mJykms~ 1 as the typical lower limit 
for our detections, the simulation gives 70% completeness 
for 400 km s -1 lines, 90% completeness for 200 km s -1 
lines, and 40% completeness for 800kms _1 lines. This 
shows that the statistical algorithm is most sensitive to 



peak intensity for lines near detection thresholds. The 
average completeness for these three widths is 67%, so 
70% is representative for a set of lines with various widths 
and 600 mJykms -1 . 

A somewhat monotonic relationship between the con- 
tinuum and line fluxes in DSFGs must exist: energy 
balance requires that galaxies with little far-IR con- 
tinuum luminosity will have weak molecular emission. 
The exact relationship is unknown, but a power-law 
fit established a representative correspondence between 
S„(CO)Av and 350 ^m flux densities S„(350) for sources 
with CO detections. This fit yielded equivalent limits of 
S„(350) = 150 mJy, 90mJy, and 50mJy for 90, 70, and 
35% completeness for the median (400 km s" 1 ) linewidth, 
as shown in Figurc[2j Given the S L/ (CO)Av-S lJ (350) dis- 
tribution, the exact form of the continuum-line relation- 
ship is unimportant over the relatively small range at low 
flux densities, and a linear fit gave essentially the same 
values. Excluding the extreme high-flux point from the 
fit flattens the S v (GO)Av~S u (350) relationship, pushing 
the completeness limits to lower continuum flux densi- 
ties. Pushing the completeness limits to higher flux den- 
sities would require a steeper relationship than could be 
supported by these data with a simple model. 

Whatever the exact form of the correspondence be- 
tween S u (CO)Av and 5^(350), all of our targets have 
^(350) falling above the 70% completeness level (prob- 
ability of detection) if their CO J = 1—0 emission falls 



Table 2. Summary of observed and derived properties for the source sample. CO data for H-ATLAS J090311. 6+003906 
(Target No. c) are from Frayer et al. (2011). Entries in italics are from tentative detections and do not have robust error 

estimates. 
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binomial distribution. Taking a lower limit of a 70% de- 
tection probability, the distribution finds 7 detections as 
most likely, with 7 ± 3 detections accounting for 99% of 
the total probability. The probability of detecting just 
two sources is 0.1%. Increasing either the detection prob- 
ability or the flux limit corresponding to a given detection 
probability reduces the probability of detecting just two 
of ten sources. 

If chance alone ruled, we should therefore have de- 
tected some 2 to 5 more weak sources. Considering the 
additional four nondetections with 350 /im flux densities 
significantly above the nominal 90% completeness limit, 
this analysis points to a strong disparity between actual 
and expected detections. We conclude that chance is 
not the only reason for nondetections, but that system- 
atic effects are also important: some of the sources could 
not be detected because their redshifts are outside of 
the Zpectrometer's band, some because they have lower 
5' ! y(CO)Aw/S' l/ (350) ratios than the detected sources, 
some because their lines are weak and broad (although 
there is no sign of these even though the eye is good at 
picking out correlated channels), or some combination. 

3.2. Redshift distribution 

Zpectrometer redshifts for the H-ATLAS sources pro- 
vide an independent sample for comparison with previ- 
ous spectroscopic redshift surveys. Comparing surveys, 
we find that the redshift distributions of the "350 /im 
peaker" galaxies with Zpectromet er detections and those 
of the radio-preselected SMGs ([Chapman et al.l 120051 ) 
that fall within the Zpectrometer's redshift range are 
strikingly similar, although the source selection and 
lines used in the redshift measurements were quite dif- 
ferent: our galaxies are 350 /im-bright targets from a 
wide-area survey that highlights lensed sources, while 
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Figure 2. Integrated CO flux Sv(CO)Av versus 350 (im flux den- 
sity 5^(350) for CO-detected and CO- undetected sources from our 
sample. CO nondetections have been assigned zero flux. Horizon- 
tal dashed lines show completeness levels for line detections de- 
rived from our 400kms — 1 linewidth simulations, with correspond- 
ing 5^(350) derived from the power- law fit shown by the smooth 
curve. All galaxies with CO nondetections have 350 (im flux den- 
sities above the 70% completeness limit for detections if a simple 
CO-continuum flux scaling law holds (see text). This implies that 
chance alone cannot explain the fraction of galaxies in the sample 
that we do not detect, but that some galaxies have redshifts out- 
side the Zpectrometer band or contain relatively little molecular 
gas. Points are labeled by target numbers given in the tables. 

within the Zpectrometer frequency range. We estimated 
the number of targets we might expect to have missed 
due to faint CO flux alone by considering the targets 
below Sj,(350) = 150 mJy, the approximate 90% com- 
pleteness level. There are 10 sources in our target list 
below this limit, of which we detect only two. If chance 
alone dominates, the detection rate will be given by the 
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the iChapman et al.l (|2005f) galaxies must be bright at 
850 £im and 20 cm radio continuum. Figure [3^ gives 
the binned distribution for the galaxies with Zpectrom- 
eter CO detections. The median of the Zpectrome- 
ter CO redshifts is z — 2A7 ± 0.11, and the mean is 
z = 2.60±0.10 (68% confidence levels by bootstrap anal- 
ysis), both well below the band center at z = 2.78. Fig- 
ure [3ji shows that the peak of the observed density func- 
tion is near z ~ 2.3, agreeing well with the peaks found 
for opti cal spectroscopic red shifts of SCUBA-selected 
sources (jChapman et al.ll2005l ). and for photometric red- 
shifts of both fferscfo d-selected sources bright at 350 /im 
(jAmblard et al.|[2010l) and 870 jum-selecte d sources iden- 
tified by LABOCA (IWardlow et all[20ll . 

Figure [3Jd gives the cumu lative distribution fun ctions 
of the Zpectrometer and the Chapm an et al.l ([20051) sam- 
ple redshifts, allowing a clean comparison free of binning 
effects. It shows that the distributions are indistinguish- 
able; a Kolmogorov-Smirnov test gives a probability of 
0.995 that the two underlying distributions are the same. 
Such good agreement must in part reflect random chance 
in a statistically small sample, but nevertheless it is clear 
that the two distributions are very similar. There is oth- 
erwise no a priori reason that t he distributions should 
be so similar: for instance, the IChapman et al.l ([20051 ) 
sample could be concentrated to lower redshifts because 
of the radio pre-selection, while the Zpectrometer detec- 
tions could highlight a population of galaxies with strong 
molecular but very little rest-frame UV line emission. 

Another similarity between the Zpectrometer sample 
and galaxies detected in millimeter-wave CO followup 
from optical redshift catalogs is the distribution of CO 
linewidths. Linewidths trace the dynamics, and thus 
to some extent masses, of the galaxies, and are not 
strongly affected by lensing. Full width at half maxi- 
mum (FWHM) linewidths in the Zpectrometer sample 
range from 210 to 1180 kms -1 , with a median width of 
400 ± lOOkrns -1 (68% confidence levels by bootstrap 
analysis). In general this is somewhat narrower than 
the w idths of mid - J line s in the millimeter - wave stud- 
ies of iGreve et all ([2005D and iTacconi et al.l (|2006l) , but 
a permutation test shows that the difference in median 
widths between the studies is significant at only the 
70% level, and could easily be from small number statis- 
tics. Combining the samples, a characteristic width of 
about 500 km s -1 is representative for DSFGs. In abso- 
lute terms, the widest 1 — lines are very broad, how- 
ever, including two with FWHM linewidths greater than 
llOOkms -1 . If tracing virialized matter, these widths 
correspond to emission from very massive galaxies or in- 
teracting massive galaxies. 

3.3. Continuum properties 

Spectroscopic redshift measu rements un ambiguously 
break the Tjj-z degeneracy (e.g.,,Blain 1999|) that renders 
dust temperature estimates uncertain for galaxies with- 
out firm redshift measurements. We fit a simple optically 
thin single-temperature dust model with dust emissiv- 
ity /? = 1.5 to the 250, 350, and 500 ^m SPIRE flux 
densities. To estimate the observed total infrared (rest- 
frame 8-1000 ^m) luminosity, we joined the far-infrared 
fit smoothly in slope to a power-l aw spectrum with form 
S v oc j/ -1 - 4 at short wavelengths (Blain et~aT1l200l . As- 
suming that the same gravitational magnification \x ap- 
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Figure 3. Comparisons of the redshift distribution of galax- 
ies from our Zpectro meter CO survey and those from the 
IChapman et ail 1 120051) optically-derived redshift survey of SMGs 
that fall within the Zpectrometer's redshift range. Panel (a) is 
the binned d istribution with Zpectrometer data in close hatch and 
IChapman et ail 1120051 ) data in coarser hatch. Panel (b) gives the 
cumulative distribution functions, which do not rely on bin width 
choices; Zpectrometer data are solid points connected by lines, 
IChapman et ail ((2005) data open triangles. The two distributions 
are strikingly similar (a Kolmogorov-Smirnov test returns a prob- 
ability of 0.995 that they are drawn from the same parent popula- 
tion) although the selection criteria are completely different. This 
suggests a common parent population (and perhaps selection ef- 
fects) for the SPIRE- and 870 /^m-selected galaxies. 

plies to emission in all far-IR wavebands, lensing should 
not affect temperature estimates, but it will scale the in- 
trinsic luminosity Lir so the observed infrared luminosity 
is ixL m . 

Table [3] contains the model results and Monte-Carlo 
error estimates for individual galaxies. Excluding the 
two z > 3.1 galaxies, the dust temperature is 34 ± 2K 
averaged over all of the detected sources, with a slight de- 
pendence on the observed 350 fim flux density (Fig. 0b) 
or, equivalently, //Ltr. Little temperature scatter is to be 
expected among the detected galaxies because the source 
selection criteria favored similar SPIRE-band SEDs. 
Temperatures near 35 K are similar t o those derived from 
other surveys of SMGs (35 ± 3 K, iKovacs et all [20061: 
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Figure 4. Panel (a) gives the observed IR luminosity versus ob- 
served 350 /xm flux density S v (ZfXi), showing that the luminosity 
scales reasonably well with flux density, even in the presence of a 
range of redshifts and dust temperatures. Panel (b) shows that 
the derived characteristic dust temperatures are nearly constant 
with iSj,(350), possibly at least in part from selection effects. The 
highest two points in both panels are from the two galaxies in our 
sample with z > 3.1. Points are labeled by target numbers given 
in the tables. 

34± 5 K. [Chapman et alJIMtOl 37 ± 1 K. IWardlow et al.l 
1201 11 1. The two highest- z galaxies in our sample have 
temperatures of about 40 K, a temperature higher than 
scatter alone can explain. Varying the dust emissivity 
parameter (3 from 1.2 to 1.7 did not change the dust 
temperatures by more than 3 K. Dust emission models 
that allow for dust emission optical depths yield temper- 
atures about 15 K higher than those from the optically 
thin model; we quote the optically thin results to facili- 
tate comparisons with previous work, most of which uses 
the same formalism. 

Figure |4ji shows that the 350 fim flux density is an ex- 
cellent proxy for the observed infrared luminosity /^Lir at 
z ~ 2—3. The two points falling above the general trend 
are the two galaxies with the highest redshifts, z > 3.1, 
in our sample. Excluding these two galaxies, we derive 
(j,L m = (5„(350 /rai)/(51mJy) +0.4) x 10 13 L o . Apart 
from the two galaxies at z > 3.1 there is no dependence 
of /xLift on redshift within errors. 

3.4. Comparison of photometric and line redshifts 

One of the goals of this project was to provide a data 
set suitable for evaluating the precision of photomet- 
ric redshifts against spectroscopic measurements. Here 



we make a first-cut comparison of techniques based on 
simple SPIRE 3-band colors and on fitting to template 
galaxy SEDs. 

With galaxies chosen as "350 /im peakers," the sim- 
plest color selection is a flux ratio that estimates the 
SED's peak wavelength. Since the 250 fim and 500 /im 
bands straddle the emission peak their ratio provides 
this estimate, reflecting some combination of tempera- 
ture and redshift. As section 13.31 shows, our sample of 
galaxies has a small range of dust temperatures, so we 
can test whether there is a simple relationship between 
the observed SED peak wavelength and redshift. Fig- 
ure [5] shows that there is a relationship between the peak 
wavelength, as given by the 250/im/500^m flux density 
ratio, and spectroscopic redshift for galaxies with CO de- 
tections. For these galaxies the 250 /im/500 /im flux den- 
sity ratio predicts redshifts within Az = 0.3 (standard 
deviation of the redshift error). Other combinations of 
continuum flux densities and ratios are less successful 
at predicting redshifts than the peak wavelength. For 
example, both the CO-detected and nondetected target 
galaxies fall along a common locus in a 500 /im/350 fim- 
350 /im/250/im color-color diagram, but the galaxies are 
well mixed in redshift along that locus. 

A more general method of estimating redshifts is to 
take SEDs of galaxies with known redshift and temper- 
ature as templates, then find which redshifted template 
best matches the observed SED. Template fits are not 
very tightly constrained by SPIRE data alone because 
the 250, 350, and 500 fim bands lie near the peak of the 
observed-frame SED, however. Flux ratios near the peak 
have little dynamic range and consequently cannot pro- 
vide strong constraints: the maximum S u (350) jS v (500) 
and 5^(350)75^(250) ratios are about 1.6. 

Figure [6] shows the redshift error Az versus redshift 
z for the detected galaxies comparing results from the 
simple 5„(250)/5,,(500) ratio and those from indepen- 
dent template fitting codes by co-authors Clements, 
Gonzalez-Nuevo, and Wardlow. Each estimates the red- 
shift by minimizing the chi-squared deviation of the data 
points compared with one or more observed SEDs of 
galaxies with known redshifts. Minimizing chi-squared 
by shifting the templates in wavelength gives the red- 
shift, while the widths of the chi-squared values ver- 
sus shift provide redshift error estimates. Comparing 
the template fitting with the linear fit between z and 
5^(250)75^(500) shown in Fig. shows that the lin- 
ear fit gives the lowest dispersion in Az for most of the 
detected galaxies. Redshifts from an SED correspond- 
ing to the Cosmic Eyelash galaxy SMM J2135— 0102 
(iGonzalez-Nuevo et al.| l2012t with SED tem plate from 
llvison et al.l 120101 and iSwinbank et al.l[20Tol ) are nearly 
as accurate as those from the linear fit. Comparison 
with models that choose from libraries of extragalactic 
SEDs shows these have generally poorer agreement than 
constrained models for most of the detected galaxies, 
with Az rs 0.5-1 and internal error estimates smaller 
than the actual redshift deviation between observation 
and model. Their errors are lower than those of the 
5^(250) / 5„(500) ratio or Eyelash fits for the two galaxies 
at z > 3.1, however. This result emphasizes the difficulty 
of fitting curves with only a few samples near the peak; 
photometric redshifts incorporating substantially longer- 
wavelength data are substantially more accurate (see, 
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Figure 5. Galaxy rcdshift z versus 250 /im/500 fim continuum 
flux density ratio (solid points; open point for a tentative detec- 
tion). Dashed horizontal lines mark the Zpcctrometer's band edges. 
CO nondetections (triangles) are plotted outside the band at an ar- 
bitrary z = 2. For galaxies with CO detections, the figure shows 
that SPIRE colors give Az ~ 0.3 for this carefully filtered homo- 
geneous sample. Points are labeled by target numbers given in the 
tables. 
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Figure 6. Photometric redshift errors versus spectroscopic red- 
shift. For the sources with CO J = 1 — redshifts, either the linear 
fit (triangles; Fig. or template fitting against the Eyelash SED 
(crosses) provides redshift estimates within Az fti 0.3 for z < 3.1. 
Dashed lines highlight grouping by method only. Template fits 
that automatically select from a library of SEDs (squares) gener- 
ally predict rcdshift to Az as 0.5-1; the grey background shows 
the typical internal uncertainty for the template fits. Models that 
choose between multiple SEDs predict redshifts outside the Zpec- 
trometer band for several of the nondetected sources, even though 
the linear and Eyelash templates place them within the band, and 
have lower errors for the z > 3 sources. As discussed in the text, 
this indicates that selection effects are important in the agreement 
between Eyelash SED and CO detections. 

e.g., ph otometric redshifts including MAMBO 1200 /im 
data in lNegrello et alJfeOlOj) . 

We emphasize that the correspondence between CO 
J = 1 - detections and SPIRE 250 Aim/500 /im color 
or Eyelash template redshifts yields no redshift informa- 
tion about the galaxies that we do not detect in the CO 
J = 1 — line. Galaxies with and without detections 
span nearly the same 250 /im/500 fim flux density ratios 
(Fig. [5]), and the success of the Eyelash template may 
simply arise from its ability to identify the T ~ 35 K 



galaxies we find in CO. As we discussed in section 13.11 
some fraction of galaxies with continuum properties sim- 
ilar to galaxies with CO detections very likely fall outside 
the Zpectrometer's band. 

3.5. Lens magnifications 

Lens magnifications (p) are needed to convert the ob- 
served fluxes to intrinsic luminosities. Determining mag- 
nification is generally done by constraining a model of 
the source-lens pair with an image of the gravitationally- 
distorted source galaxy. Detailed lens models are not yet 
available for most of our sources. We do have velocity- 
resolved spectra, however, and can take an alternative 
approach of comparing observed and intrinsic luminosi- 
ties established by an empirical luminosity-linewidth re- 
lationship similar in spirit to the Tully-Fisher relation 
(Tull y fc Fisherl 119771 ) . Without correction for galaxy 
inclination or dispersion in intrinsic galaxy properties, 
such a relationship cannot be exact, but it can provide 
approximate estimates of lensing magnifications. 

Since the galaxies with Zpectrometer CO detections 
seem to be quite typical of the general SMC population, 
as we discussed above, we assume that only lens magnifi- 
cations modify the observed typical luminosity distribu- 
tions of galaxi es detected in CO J = 1 — 0. Following 
the method of Bothwcll et al. (2012), as outlined below, 
we find an empirical intrinsic integrated line luminosity- 
linewidth relationship L' = a (A^fwhm) 6 from galaxies 
with published CO J — 1 — intensities, widths, and 
magnifications. This relationship forms the basis to solve 
for the unknown magnification /i that scales the true line 
luminosity L' to an apparent luminosity ^apparent = pL', 
or 



H = 



apparent 



L' 



apparent 



(i) 



a (A« FWH m) 

Figure shows £coj=i-o for 15 z = 2 ~ 4 SMGs 
from the literature versus linewidth, with corrections 
for lens magnification when needed (CO data from 
Harris et all [20101: ICarilli et all IMA llvison et aT1l20Tl1: 
Riechers et all 1 20 11 at magnifications as needed from 
Smail et all 120021 ) . Points for these galaxies fall close 
to a power law fit with relatively small scatter. Inserting 
the fit parameters into equation ([T]), we obtain 



/i = 3.5 



Li 



apparent 



10 11 Kkms' 1 pc 2 



400kms~ J 

AfFWHM 



1.7 



(2) 



with the unit s for line luminosity L' as defined in 
ISolomon et all (|1992l ). Apparent line luminosities for 
the H- ATLAS galaxies with Zpectrometer CO measure- 
ments, uncorrected for lens magnification, are also plot- 
ted in Figure [7] . As expected for a flux-limited sam- 
ple, these galaxies have a narrower range of luminosities 
than the comparison sample from the literature. Galax- 
ies with narrower line widths, which will generally have 
lower luminosities, require more magnification to reach 
the observational detection threshold. 

Table |3] lists the lens magnifications derived from equa- 
tion The empirical magnification distribution is ap- 
proximately uniform over its range of 3 to 20, with both 
median and mean equal to 10. Simple uncertainty es- 
timates for the scaling term and power law index are 
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Figure 7. Square points indicate line luminosities L'^q ver- 
sus CO J = 1 — FWHM linewidths for 15 sources with red- 
shifts, flux es, and magnificat i ons as needed from the literature (CO 
data fr om Harris ct al. 2010; ICarilli et aT|[20TOl; llvison et ani20TTl; 
I Hiechers et al. 20 llal ; magnifications as needed from ISmail et alj 
2002). The solid line shows a power-law fit to these data, represent- 
ing an estimate of the intrinsic line luminosity versus linewidth for 
typical galaxies. Round points show the observed line luminosities 
versus widths for the Zpectrometer sample (values from Table [SJ, 
without correction for lens magnifications. The literature sample 
has a steep power-law relation with modest scatter, while the Zpec- 
trometer sample shows a much flatter trend. Intrinsically weaker 
sources require higher magnification to rise above the observational 
detection threshold, which is approximately constant in flux. We 
scaled the observed line luminosity for each of the Zpectrometer 
sources by the prediction from the power-law estimate of its likely 
magnification fi, as given in Table [3] This method yielded a range 
of fi 3-20 and a median of fi = 10 for the Zpectrometer sample. 

not available because the two are not independent, but 
a bootstrap analysis provides ^68% confidence bounds 
on /i for each galaxy in Table El Fractional uncertainties 
range from about 30% at a minimum near 500 km s -1 , 
climbing to about 50% at 300 and 1200kms _1 , and 
about 100% at 200kms _1 . The divergence toward low 
A^fwhm is also a reminder that the (unknown) incli- 
nation corrections become large for galaxies with ob- 
served low linewidths. It is possible that some of the 
sources with broad CO lines and low inferred amplifica- 
tions from equation ^ are intrinsically hyper-luminous 
systems which are not significantly magnified at all; the 
magnification uncertainty encompasses that possibility. 

Magnifications from the luminosity-linewidth relation 
([2]) compare well with values from the two galaxies in 
the H- ATLAS survey that have both Zpectrometer de- 
tections and preliminary lens models from image recon- 
struction with lens and source galaxies at known red- 
shifts. For J0903 11. 6+003906 (ID.81; CO data from 
iFrayer et al.1 120111 ). the CO linewidth-luminosity rela- 
tionship in equation ([2]) predicts /i = 14 ± 4, while 
INegrello et aiT(|2010D quote fx ~ 19 (range 18-31) from 
their modeling. For J091305.0-0 05343 (ID. 130) equa- 
tion (0) predicts /i = 10 ± 4 while INegrello et al.1 ()2010f ) 
quote n ~ 6 (range 5-7). 

Table [3] provides estimates of the intrinsic luminosities 
and gas masses corrected by the source-specific magnifi- 
cation factors, and also by the median /i = 10. Most of 
the galaxies appear to have intrinsic luminosities typical 
of ULIRGs, although three are candidate hyper-LIRGs, 
with luminosities greater than 10 13 Lq. Gas mass esti- 



mates use the same CO lumi nosit y-to-mass conversion 
factor as in IGreve et al.1 (f2005h and lTacconi et~aTI §006) 
for ease of comparison, a = 0.8 M Q (K kms _1 pc 2 ) 
(|Downes & Solomon! 1 1 9981 ) . Although absolute mass es- 
timates carry double uncertainties from both conversion 
and magnification factors, they are similar to estimates 
for SMGs selected from 850 /im observations that use the 
same conversion factor. 

4. DISCUSSION & CONCLUSIONS 

Demonstrating an ability to rapidly establish spectro- 
scopic redshifts and molecular line parameters for lumi- 
nous high-redshift galaxies is a major step forward in 
exploring the properties of DSFGs. In a single observing 
season we have added a dozen galaxies with CO J = 1 — 
line detections to the pool of sources for study among the 
bright "350 fim peakers" from the Herschel- ATLAS sur- 
vey. 

We established accurate redshifts for sufficient galax- 
ies to compare with an independent survey of SMGs with 
observed optical redshifts. That survey started with sub- 
millimeter source continuum detections in ~ 15 arcsec- 
ond beams at 850 /im, associated those with 20 cm ra- 
dio continuum sources to obtain positions to ~ 1 arc- 
second, and then used optical spectroscopy to deter- 
mine redshifts (see e.g.. iChapman etlrtl 120051 ). Subse- 
quent observations with a millimeter-wave interferome- 
ter used the optical redshifts as starting points to deter- 
mine molecular redshi f ts from the star-forming gas itself 
(e^ iNeri et~a^[200l IGreve et al.l [20051: ITacconi et all 
2006), conclusively connecting the 850 /im sources with 
massive molecular clouds in the target galaxies in about 
half of the attempts. The Zpectrometer CO redshifts, 
in contrast, started with bright 350 /im sources with 
submillimeter photometry indicating rough photometric 
redshifts, but then made a single step from continuum 
to molecular line observations with this wideband spec- 
trometer. We find nearly identical redshift distributions 
for the two surveys, indicating that each sample is rep- 
resentative of a general DSFG population. Coupled with 
mm- wave molecular s pectroscopy of many of the brighter 
850 /im source s (e.g ., (Neri et al.ll200l IGreve et al.1 12005: 
ITacconi et al"1l2006f ). the agreement in distributions be- 
tween independent samples lays to rest some of the con- 
cern that SMG redshift distributions from optical spec- 
troscopy may be influenced by misidentifications between 
the submillimeter sources and counterpart optical galax- 
ies. 

Dust temperatures and linewidths of the galax- 
ies with CO 1-0 detections appear to be similar to 
SMGs. With spectroscopic redshifts, we are able 
to break the temperature-redshift degeneracy to find 
a characteristic temperature of 34 K (in an optically 
thin formal ism) , similar to those typically reported 
for SMGs (IKovacs et al] l200l IChapman etldl 12011 
IWardlow et al.ll2011l) . Galaxies with redshifts z > 3 may 
be somewhat warmer than those at lower redshift. The 
CO linewidths have very similar distributions for these 
H-ATLAS CO J = 1 - detections and for SMGs ob- 
served in mid-J lines. Both of these comparisons add 
further to an argument that the bulk of samples selected 
by bright 850 /Ltm and SPIRE emission draw from the 
same parent population. 
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Table 3 

List of galaxies with observed CO line luminosities and estimates for lens magnifications fi and corrected 

(intrinsic) dust luminosities and gas masses. The median value of fi = 10. Mass calculations are from 
M = a L', with a = 0.8 M (K km s" 1 pc 2 ) -1 IIDownes fc So lomon 199D to ease comparisons using this 
value as a standard. Masses from J = 1 — intensities are likely to provide larger masses by a factor between 

1.5 and 2 comp ared with those derived from millimeter-wave obser vations (discussion a nd references in 
Harris ct al. 2010). Errors do not include uncertainties in a (see, e.g.. ITacconi et al.H2008l . for a review). For 
comparison with other work that assumes /i ~ 10, the last two columns contain luminosity and mass 
estimates with fi = 10, with errors incorporating a nominal Afi//i = 0.5. 



H-ATLAS Name 


^bs.co I 1010 K 


l-i 


Ljr/h 


M gas //i 


L IR /W 


M gas /10 


Target 




km s 1 pc 2 ] 




[10 12 L e ] 


[10 10 M Q ] 


[10 12 L Q ] 


[10 10 M ] 


No. 


J084933.4+021443 


29.2 ± 10.4 


2 + 1 


28.8 + 19.6 


11.7 + 8.0 


5.8 + 3.5 


2.3 + 1.4 


2 


J090302. 9-014127 


26.1 + 3.3 


18 + 8 


3.8 + 1.9 


1.2 + 0.6 


6.9 + 3.5 


2.1 ± 1.1 


(i 


J090311. 6+003906 


46.9 + 3.8 


14 + 4 


4.1 ± 1.2 


2.7 + 0.8 


5.8 + 2.9 


3.8 + 1.9 


c 


J091305.0-005343 


24.9 + 3.8 


10 + 4 


3.1 ± 1.2 


2.0 + 0.8 


3.1 ± 1.6 


2.0 + 1.0 


b 


J091840. 8+023047 


32.9 + 8.3 


5 + 1 


8.0 + 2.9 


5.3 + 1.9 


4.0 + 2.2 


2.6 + 1.5 


5 


J113243. 1-005108 


20.7 + 6.0 


8 + 3 


3.0 + 1.3 


2.1 + 0.9 


2.4 + 1.4 


1.7 + 1.0 


10 


J113526. 3-014605 


15.4 + 3.3 


17 + 11 


7.1 + 1.8 


0.7 + 0.5 


12.1 + 6.6 


1.2 + 0.7 


9 


J114637.9-001132 


46.7 + 7.4 


7 + 2 


22.4 + 6.9 


5.3 ± 1.6 


15.7 + 8.2 


3.7 + 2.0 


13 


J115820. 2-013753 


17.6 + 2.8 


13 + 7 


1.9 + 1.0 


1.1 + 0.6 


2.5 ± 1.3 


1.4 + 0.7 


15 


J133649. 9+291801 


22.2 + 2.9 


23 + 15 


2.5 ± 1.6 


0.8 + 0.5 


5.7 + 2.9 


1.8 + 0.9 


20 


J134429. 4+303036 


71.2 ± 10.0 


4 + 2 


25.2 ± 14.2 


14.2 ± 8.0 


10.1 ± 5.2 


5.7 + 3.0 


19 


J141351. 9-000026 


43.4 + 5.0 


10 + 3 


4.9 + 1.3 


3.5 + 0.9 


4.9 + 2.5 


3.5 + 1.8 


21 



Most, if not all, of the galaxies we observed appear 
to be gravitationally lensed. We estimate lens mag- 
nifications with an empirical CO linewidth-luminosity 
relationship made possible with velocity-resolved spec- 
troscopy. We verified magnification agreement within 
factors close to two for two galaxies with published 
lens modeling. We find magnifications from 3 to 20, 
a range matc hing those for lenses with detailed mod- 
els, m ~ 5-30 JSwinbank et al.llMot INegrello et aT1 l2TTifl 
IGavazzi et al.|[2*011l ; ILestrade et al.ll2011[ ). In the absence 
of other information, a typical value of /i = 10 is likely to 
be correct within a factor of two for the bright H-ATLAS 
galaxies. 

With even moderately accurate lens magnification es- 
timates, we can make reasonable estimates of intrinsic 
luminosities and gas masses of the galaxies in our sam- 
ple (Table [3]) . Most have luminosities characteristic of 
ULIRGs, but there are three targets that could be hyper- 
LIRGs with luminosities > 1O 13 L . The most luminous 
of these could well be multiple galaxies within the beam, 
as their linewidths are broad. Even accounting for lens- 
ing magnification, however, the bright H-ATLAS galax- 
ies seem to be drawn from the high-luminosity, high-mass 
end of the DSFG distribution. 

The overall agreement in properties between the Zpec- 
trometer and 850-1200 /mi-selected SMG samples indi- 
cates that the two populations are drawn from the same 
distribution. In principle, the additional selection in- 
troduced by lens magnification could have emphasized 
different properties from the SMG population, but it 
does not appear to have done so. Lensed sources will 
be valuable in understanding the properties of the un- 
derlying DSFG population, including the SMGs: for ob- 
servations, magnification by even a factor of a few trans- 
lates to an order of magnitude reduction in integration 
time. Detailed studies of the galaxies' interstellar media 
in multiple molecular and atomic lines becomes feasibl e 
in the bright lensed sample (e.g.. iDanielson et al.ll201l[ ). 
with magnification al so increasing spatial r esolution in 
the source plane (e.g.. ISwinbank et al.H2010l ). 



Overall, the detection fraction of about half of the tar- 
gets compares well with other programs that observed 
molecular lines from high redshift galaxies. A conserva- 
tive extrapolation of a CO-to-continuum intensity rela- 
tionship indicates that the nondctcction rate is substan- 
tially higher than expected from chance alone, evidence 
that some of the simple photometric redshifts are incor- 
rect. Our selection is from SPIRE photometry identifi- 
cation of "350 /mi peakers." Comparison of the SPIRE 
data and CO redshifts indicates that the continuum peak 
wavelength, deduced from the 250 /im/500/im flux den- 
sity ratio, predicts the redshift of the galaxies from our 
carefully-selected sample with Az ~ 0.3. The SPIRE 
colors do not necessarily provide accurate estimates for 
galaxies we do not detect, however. A simple explana- 
tion for some of the nondetections is emission from dust 
warmer than the ~ 35 K characteristic for most of the 
galaxies with CO detections. The redshift-tcmpcraturc 
degeneracy for the observed peak emissio n wavelength, 

^peak,obs C*- 

T/(l + z) (e.g., Hihil 1999) causes prob- 
lems for photometric selection in the SPIRE bands alone: 
warmer galaxies with redshift above the Zpectrometer's 
upper limit will appear in the target list. It is very likely 
that this is the reason for the nondetection of the target 
with third-brightest 350 /xm flux density in this study. 

Reasons other than warm galaxies or CO intensities 
that fall below our observational limit can also contribute 
to nondetections. Emission from multiple galaxies within 
the beam at one or more wavelengths can contribute flux 
to distort the apparent SED, leading to erroneous red- 
shifts from single-component models. This is likely to 
affect weaker sources preferentially, both because smaller 
absolute amounts of contaminating flux are increasingly 
important and because the density of contaminating 
sources increases with decreasing flux. Multiple distinct 
temperature components within a source, for instance a 
central AGN-heated dust region surrounded by less ex- 
cited material, will also distort the SED. Depending upon 
the temperature structure within the source and the size 
and position of the source relative to any lens caustic, 
differential amplification could distort observed SEDs. 
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SEDs in all these cases will be poorly represented by 
the simple single-component model we have used in tar- 
get selection, and photometric redshifts based on SPIRE 
data alone will fail to some extent. 

The observed redshift distribution and the dust SED 
properties hint at some change in DSFG properties near 
z ~ 3. The redshift cumula tive distribution f unctio ns for 
both the Zpectrometer and lChapman et~a l. (2005) sam- 
ples show a slope break at z ~ 2.9 although they are 
selected in different ways. This indicates that DSFGs 
become less common, or are selected less efficiently by 
continuum properties, above that redshift. An increase 
in the number of warm galaxies with redshift would both 
explain the tentative trend in temperature we see and 
explain the lower CO detection rate. Temperature alone 
may not suffice to explain the break, and explanations 
involving more compl ex SEDs with increas ing redshift 
are attractive as well. IWardlow et al.l (|2011[ ) use exten- 
sive optical and infrared photometry on a sample of 74 
870 ^mi-selected SMGs, also finding a precipitous drop in 
submillimeter source counts above z p hot ~ 2.8 in spite of 
their sensitivity to emission from warm dust. Increased 
sample sizes are necessary to establish whether there is 
a real change in DSFG properties near z ~ 3 or not. 

The CO data presented here are among the earliest 
spectroscopic explorations of Herschel-selected galaxies. 
Continuing observations of the CO J = 1 — line by the 
GBT and Jansky-VLA will provide a more diverse sam- 
ple of galaxies: some from broadened continuum selec- 
tion criteria, and some with deeper integrations to reach 
weaker line emission. Wideband spectrometers now com- 
ing on line for the 3 mm band at the LMT and IRAM 
30 m will also add to the number of galaxies with known 
molecular redshifts. An expanded sample will strengthen 
our understanding of the relationship between molecular 
and continuum emission in DSFGs. 
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